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ABSTRACT 

We investigate the calibration and uncertainties of black hole mass estimates based on the single-epoch 
(SE) method, using homogeneous and high-quality multi-epoch spectra obtained by the Lick Active Galactic 
Nucleus (AGN) Monitoring Project for 9 local Seyfert 1 galaxies with black hole masses < 10*^ Mq. By 
decomposing the spectra into their AGN and stellar components, we study the variability of the single-epoch 
Hf3 line width (full width at half-maximum intensity, FWHMh/3 ; or dispersion, cth/?) and of the AGN continuum 
luminosity at 5100 A (Lsioo)- From the distribution of the "virial products" (oc FWHMh/3^ ^sioo cth/S^ 
Ljjqq) measured from SE spectra, we estimate the uncertainty due to the combined variability as ~ 0.05 dex 
(12%). This is subdominant with respect to the total uncertainty in SE mass estimates, which is dominated by 
uncertainties in the size-luminosity relation and virial coefficient, and is estimated to be ~ 0.46 dex (factor of 
~ 3). By comparing the H/3 line profile of the SE, mean, and root-mean-square (rms) spectra, we find that the 
H/3 line is broader in the mean (and SE) spectra than in the rms spectra by ^0.1 dex (25%) for our sample 
with FWHMh^ < 3000 km s"'. This result is at variance with larger mass black holes where the difference is 
typically found to be much less than 0.1 dex. To correct for this systematic difference of the H/3 line profile, 
we introduce a line-width dependent virial factor, resulting in a recalibration of SE black hole mass estimators 
for low-mass AGNs. 

Subject headings: galaxies; active - galaxies: nuclei - galaxies: Seyfert 
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L INTRODUCTION 

■ Supermassive black holes (BHs) are believed to play a key 
^ role in galaxy evolution. Evidence for this connection comes 
• ^ from the tight correlations observed in the local universe be- 
tween BH masses and the global properties of their host galax- 
^ ies (Magon-ian et al. 1998; FeiTai-ese & Memtt 2000; Geb- 
. 5^ hai-dt et al. 2000; Gultekin et al. 2009; Bentz et al. 2009a; 
Woo et al. 2010). Establishing the cosmic evolution of these 
correlations is a powerful way to understand the feedback 
mechanisms connecting BHs and galaxies (e.g., Kauffmann & 
Haehnelt 2000; Robertson et al. 2006; Hopkins et al. 2006). 
Recent observational studies have found that these correla- 
tions may evolve over cosmic time, in the sense that BHs of 
a given mass appeared to live in smaller galaxies in the past 
(e.g.. Woo et al. 2006; Peng et al. 2006; Treu et al. 2007; Woo 
et al. 2008; Merloni et al. 2010; Decarli et al. 2010; Bennert 
et al. 2010). 

In order to investigate the nature of BH-galaxy coevolu- 
tion, as well as virtually all aspects of active galactic nu- 
cleus (AGN) physics (e.g., Woo & Urry 2002; Kollmeier et 
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al. 2006; Davis et al. 2007), BH masses must be accurately 
determined at large distances. Dynamical methods based on 
high angular resolution kinematics of stars and gas are the 
most common approach to measuring masses of quiescent 
BHs (e.g., Kormendy & Gebhardt 2001; Ferrarese & Ford 
2005). However, owing to the parsec-size scale of the sphere 
of the influence of typical BHs, they are limited to galaxies 
within a distance of few tens of Mpc with current technology. 

In the case of BHs powering an AGN, the presence of a 
variable broad-line region (BLR) provides an alternative way 
that is in principle applicable to much larger distances. The 
geometry and kinematics of the BLR gas can be mapped in the 
time domain using the so-called reverberation (or echo) map- 
ping technique (Blandford & McKee 1982; Peterson 1993). 
In turn, these quantities can be converted into BH mass esti- 
mates under appropriate assumptions about the dynamics of 
the system (Peterson 1993; Pancoast, Brewer, & Treu 2011). 
Estimators of the form Mbh oc Rblr V^, where Rblr and V 
are (respectively) size and velocity estimators of the BLR, 
are often referred to as "virial" mass estimators. However, 
due to the observational challenges of reverberation mapping 
campaigns, fewer than 50 BH masses have been measured to 
date using this technique (Wandel, Peterson, & MaUcan 1999; 
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Kaspi et al. 2000; Peterson et al. 2004; Bentz et al. 2009c; 
Denney et al 2009; Barth et al. 2011). 

In light of the scientific importance of determining BH 
masses, it is critical to develop alternative BH mass estima- 
tors that are observationally less demanding. A popular BH 
mass estimator, based on the results of reverberation mapping 
studies, is the so-called single-epoch (SE) method. It exploits 
the empirical correlation between the size of the BLR and the 
AGN continuum luminosity (/?blr oc L", with a k, 0.5), as 
expected from the photoionization model predictions (Wan- 
del, Peterson, & Malkan 1999; Kaspi et al. 2000, 2005; Bentz 
et al. 2006, 2009b), to bypass the expense of a monitoring 
campaign. Thus, the AGN luminosity is used as a proxy for 
the BLR size and, in combination with the square of a ve- 
locity estimate from a broad line, to estimate BH masses from 
single spectroscopic observations. Typically, SE mass estima- 
tors are based on optical/ultraviolet lines (e.g., H/3 or Mg II) 
and optical/ultraviolet continuum luminosity (e.g., at 5100 A 
or 3000 A). A summary and cross-calibration of commonly 
adopted recipes is given by McGill et al. (2008). 

Due to its convenience, the SE method has been widely ap- 
plied from the study of BH demographics (e.g., Shen et al. 
2008; Fine et al. 2008) to the characterization of galaxy- 
AGN scaling relations at low and high redshift (e.g., Treu et 
al. 2004; Barth et al. 2005; Greene & Ho 2006; Woo et al. 
2006; Bennert et al. 2010, Bennert et al. 2011a). For this rea- 
son it is of paramount importance to quantify, understand, and 
(possibly) correct for random and systematic uncertainties in 
the method. In addition to the random and systematic errors, 
selection bias can play a role in studying statistical properties 
of AGN samples selected from a flux-limited survey since BH 
mass from SE data is proportional to the AGN continuum lu- 
minosity at 5100A (L5100) (e.g., Lauer et al. 2007; Treu et 
al. 2007; Shen & Kelly 2010). Naturally, the strength of the 
selection bias depends on the uncertainty of the SE mass es- 
timates, providing another compelling reason to quantify it 
accurately. 

The largest uncertainty comes from the unknown "virial" 
factor /, connecting the observable size and velocity to the 
actual BH mass, Mbh = /^blr /G, where G is the grav- 
itational constant. In general, / cannot be determined for 
individual sources due to limited spatial information except 
a few cases (Davies et al. 2006, Onken et al. 2007; Hicks 
& Malkan 2008; see, however. Brewer et al. 2011 and ref- 
erences therein). Therefore, an average virial factor is typi- 
cally applied. This average is determined by forcing active 
and quiescent galaxies to obey the same BH mass-galaxy ve- 
locity dispersion (Mbh-o"*) relation (Onken et al. 2004; Woo 
et al. 2010), even though the virial factor of individual AGNs 
may be different from the mean value. Thus, using an av- 
erage virial factor introduces an uncertainty in the SE mass 
estimates. It is not known precisely how large the uncertainty 
of the virial factor is (see Collin et al. 2006), and whether this 
uncertainty is stochastic (random) or has a systematic compo- 
nent that can be reduced using additional observables. An up- 
per limit to the uncertainty is derived from the intrinsic scatter 
of the AGN Mbh -cr* relation (0.43 dex; Woo et al. 2010), as- 
suming that the samples used to calibrate / are representative 
of the class of broad-line AGNs targeted for the SE study. 

A second source of uncertainty is the variability of AGNs: 
line width and continuum luminosity will vary as a function 
of time, while the BH mass is not expected to change sig- 
nificantly over time scales of order a few years. Thus, AGN 



variability introduces an uncertainty in the SE mass estimates, 
which is believed to be stochastic in nature. Previous studies 
based on multi-epoch spectra reported that the random error 
due to the variability is ^15-25% (e.g.. Woo et al. 2007; Den- 
ney et al. 2009). 

A third source of error is the intrinsic scatter in the size- 
luminosity relation used to infer the size of the BLR. Recent 
studies, based on reverberation mapping results and Hubble 
Space Telescope (HST) imaging analysis, report ^40% scatter 
in the size-luminosity relation (Bentz et al. 2009b). 

A fourth source of uncertainty in SE mass estimates is 
due to differences in the BLR line profile as measured from 
SE spectra and those measured from root-mean-square (rms) 
spectra. In reverberation mapping studies, BH mass determi- 
nations rely on the line width measured from the rms spec- 
tra, which reflect the varying part of the line profile. In con- 
trast, for SE mass determinations line widths are measured 
from single spectra since such equivalent measurements as 
in the rms spectra are not available. Thus, it is necessary 
to investigate and quantify the line-width difference between 
SE and rms spectra. Previous studies based on multi-epoch 
data showed that the H/3 line widths in the mean spectra are 
broader than those in the rms spectra (e.g., Sergeev et al. 
1999; Shapovalova et al. 2004; ColHn et al. 2006; Denney 
et al. 2009). The difference is presumably due to the different 
kinematics of the gas responding over various time scales, in- 
dicating that a different normalization is required in order to 
consistently estimate virial masses based on the SE method. 

In this work, we focus on the uncertainties of SE mass es- 
timators due to the variability, and those due to differences 
in line profiles. By comparing measurements from single- 
epoch, mean, and rms spectra using the high-quality multi- 
epoch spectra of 9 local Seyfert galaxies in the relatively un- 
explored regime of low-mass BHs from the Lick AGN mon- 
itoring project (Bentz et al. 2009c), we provide new quantita- 
tive estimates for the uncertainties and recipes to correct for 
them. The paper is organized as follows. In ^ we describe 
the observations and data reduction. Section [3] discusses the 
measurement method for SE spectra as well as for mean and 
rms spectra. In § H] we present the main results including a 
test of the virial assumption, a quantification of the random 
errors due to AGN variability, and the systematic differences 
in line width between SE and rms spectra. We also present a 
recalibration of standard recipes that corrects for the system- 
atic differences. We conclude and summarize our findings in 
§ ID Throughout the paper we adopt the following cosmolog- 
ical parameters to calculate distances: Hq = 70 km s~' Mpc~', 
Qn, = 0.30, and Qa = 0.70. 

2. OBSERVATIONS AND DATA REDUCTION 

We use the homogeneous and high-quality multi-epoch 
spectra from t he Lick AGN monitoring project (LAMP; 
iBentz et aP '2009c). which was designed to measure the re- 
verberation time scales of 13 local Seyfert 1 galaxies. Here, 
we briefly summarize the observations and data reduction. 

The LAMP campaign was carried out using the Kast spec- 
trograph at the 3-m Shane telescope at the Lick Observatory in 
Spring 2008. Among 13 Seyfert 1 galaxies, we selected 9 ob- 
jects for which the H/? line variability was sufficiently large to 
measure the reverberation time lag (Bentz et al. 2009c). Dur- 
ing the LAMP campaign, each object was observed multiple 
times (43 to 5 1 epochs with an average of 47), enabling us to 
construct high-quality mean and rm.s. spectra. 

After performing standard spectroscopic reductions using 
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one-dimensional spectra were extracted with an aper- 
ture window of 13 pixels (10"1). Flux calibrations utilized 
nightly spectra of spectrophotometric standard stars. As de- 
scribed by Bentz et al. (2009c), the spectral rescaUng was 
performed using the algorithm of van Groningen & Wanders 
JT992) in order to mitigate the effects of slit loss, variable see- 
ing, and transparency. By rescaling, shifting, and smoothing 
each spectrum, the algorithm minimizes the difference of flux 
of the [O III] lines between each spectrum and a reference 
spectrum created from the mean of individual spectra. The 
quality of each individual spectrum is sufficient to perform 
SE measurements (average signal-to-noise ratio S/N w 66 per 
pixel at rest-frame 5100 A). 

3. MEASUREMENTS 

Two quantities, the line width and the continuum luminos- 
ity, are required to determine Mbh using single spectra. Uni- 
form and consistent analysis is crucial for investigating sys- 
tematic uncertainties and minimizing additional errors. In 
this section, we present the multi-component spectral fitting 
process and describe the measurements using single-epoch, 
mean, and rms spectra. 

3.1. Multi-Component Fitting 
To measure the line width of H/3 and the continuum lumi- 
nosity at 5 100 A, we follow the procedure given bv lWoo et al.l 
((2006) and McGifl et al. (2008), but with significant modifi- 
cations as described below (cf., McLure & Dunlop 2004; Di- 
etrich et al. 2005; Denney et al. 2009). The multi-component 
fitting processes were carried out in a simultaneous and au- 
tomated fashion, using the nonlinear Levenberg-Marquardt 
least-squares fitting routine mpf it (Markwardt 2009) in IDL. 

First, all single spectra were converted to the rest frame. 
Second, we modelled the observed continuum with three 
components: the featureless AGN continuum, the Fe II emis- 
sion blends, and the host-galaxy starlight, using respec- 
tively a single power-law continuum, an Fe II template from 
Boroson & Green (1992), and a host-galaxy template from 
Bruzual & Chariot (2003). A simple stellar population syn- 
thesis model with solar metallicity and age of 11 Gyr was 
found to reproduce the observed stellar lines reasonably well 
(see Figure 1). The Fe II emission blends and the host-galaxy 
template were convolved with appropriate Gaussian velocities 
to reproduce kinetic and instrumental broadening during the 
fitting process as described below. The best continuum mod- 
els were determined based on the statistic in the regions 
4430-4600 A and 5080-5550 A where Fe II emission domi- 
nates. The three components were varied simultaneously with 
six free parameters: the normalization and the slope of the 
power-law continuum, the strength and the broadening veloc- 
ity of the Fe II, and the line strength and the velocity disper- 
sion of the host-galaxy templates. We masked out the typical 
weak AGN narrow emission lines (e.g.. He I A4471, [Fe VII] 
A5160, [N I] A5201 , [Ca V] A5310; Vanden Berk et al. 2001) 
during the fitting process. The best-fit continuum models (the 
power-law component + the Fe II template + the host-galaxy 
template) were subtracted from each spectrum, leaving the 
broad and narrow AGN emission lines. 

IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation (NSF). 



TABLE 1 

H/3 Integration Ranges, and H/3 narrow ratios 





(A) 


j\^Pnc)'J\1}~' mi A^uu/j 


Arp 151 


4790^980 


0.18 


NGC4748 


4790^920 


0.12 


Mrkl31() 


4800^920 


0.13 


Mrk202 


4810^920 


0.35 


NGC4253 (Mrk766) 


4790^930 


0.13 


NGC6814 


4760^950 


0.03 


SBS 1116+583A 


4795^940 


0.12 


Mrk 142 


4790^910 


0.37 


NGC5548 


4705-5040 


0.10 



Note. — All values in the table are given in the rest frame. 



Third, we subtracted the narrow lines around the H/3 region 
before fitting the broad component. We first made a template 
for the H/3 narrow-line profile by fitting a tenth-order Gauss- 
Hermite series (cf., van der Marel & Franx 1993) model to the 
[O III] A5007 line. We then subtracted the [O III] A4959 line 
by blueshifting and scaling the flux of the template by 1/3. 
The H/3 narrow line was also subtracted by scaling the [O III] 
A5007 Une. The ratios of the narrow H/3 to [O III] A5007 were 
determined from the minimization in the mean spectra and 
then forced to be the same for all SE spectra of each object. 
Applied scaling ratios for the H/3 narrow component range 
from 0.03 to 0.37 (Tableffl). 

Lastly, we modelled the broad component of the H/3 line 
using a sixth-order Gauss-Hermite series. We also used a two- 
component Gaussian model to describe the broad and narrow 
components of the He II A4686 emission line whenever it af- 
fected the blue wing of the H/3 profile. Figure [U shows the 
fitting results for the mean spectra. 

3.2. Single-Epoch Spectra 

We performed the multi-component fitting procedure using 
individual SE spectra, and measured the line width and con- 
tinuum luminosity for each epoch. The vast majority of SE 
spectra have sufficiently high quality to perform the analy- 
sis (S/N « 66 per pixel at rest-frame 5100 A). However, a 
small fraction of spectra have significantly lower S/N owing 
to bad weather during the LAMP monitoring campaign. In ad- 
dition, there are a few epochs with artificial signatures, such as 
bad pixels, abnormal curvature, or fluctuations in the reduced 
spectra. Those SE spectra were discarded to avoid possible bi- 
ases due to much larger measurement errors (see Fig. |2]). On 
average, four bad epochs out of 47 nights were removed for 
each object, except for SBS 1 1 16, for which 1 1 epochs were 
eliminated because of a defect between the H/3 and [O III] 
A4959 lines due to bad pixels in the detector. 

3.2.1. Emission-Line Width 

We measured the full width at half-maximum intensity 
(FWHMh^) a nd the dispersion (ct ha) (the second moment 
of line profile; iPeterson et al.l 120041) of the broad component 
of H/3 directly from the data as well as from the fits to the 
continuum-subtracted spectra. Line-width measurements are 
corrected for the instrumental resolution in a standard way 
(Barth et al. 2002; Woo et al. 2004; Bentz et al. 2009c), by 
subtracting in quadrature the instrumental resolution (Table 
1 1 of Bentz et al. 2009c) from the measured line width. 

By comparing line widths measured from Gauss-Hermite 
series fits with those directly measured from the data, we 
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Rest Wavelength (A) 

Fig . 1 . — Multi-component spectral fitting in the mean spectra. The mean spectra of all 9 Seyfert galaxies are presented along with multi-component models. In 
each panel, observed spectra (black) and the continuum+Fe IH-stellar best-fit model (magenta) are shown in the upper part, and the best-fit power-law continuum 
(green), stellar spectrum (yellow), and Fe II template (violet) models are presented in the middle part. Three narrow lines [H/3, [O III] AA4959,5007 (blue)], 
broad H/3 (red), and the broad and narrow He II A4686 components (brown) are presented in the bottom part. The residuals (black), representing the difference 
between the observed spectra and the sum of all model components, are arbitrarily shifted downward for clarity. 



found less than a 3% systematic difference (with consider- 
able rms scatter of ^5%) as expected given the high S/N 
of individual spectra. The small systematic trend between 
FWHMh^ and <Jufs shows opposite directions. In the case 
of FWHIVIh^, the measurements from the fit were 2.6 ± 0.2% 
larger than those from the data while measurements from 
the fit were 1 .9 ± 0. 1 % smaller than those from the data, show- 
ing a trend consistent with that reported by Denney et al. 
(2009). For consistency with other studies on the reverber- 
ation and single-epoch masses, we focus on the line-width 
measurements from the fits in the rest of the paper unless ex- 
plicitly noted. 

3.2.2. Continuum Luminosity 

We measured the monochromatic continuum luminosity at 
5100 A from the observed spectra at each epoch by calculat- 
ing the average flux in the rest-frame 5080-5120 A region. 
The luminosity at 5100 A (total luminosity, Lsioo.t) is strongly 
contaminated by the host-galaxy starlight when the AGN lu- 
minosity is comparable to or smaller than the host-galaxy stel- 
lar luminosity as in the Seyfert galaxies in our sample. 



TABLE 2 

Mean continuum luminosities and host-galaxy contributions 



Object ALA(tot) ALa(AGN) ALA(star) 

(Iff^^ergs-') (10''2 ergs-') (lO^ergs"') 



(1) (2) (3) (4) (5) 



Arpl51 6.11 2.76 3.42 0.56 

NGC4748 10.05 5.55 4.69 0.47 

Mrkl310 8.10 3.79 4.32 0.53 

Mrk202 8.78 5.02 3.75 0.43 

NGC4253 8.87 6.41 2.47 0.28 

NGC6814 1.93 0.83 1.15 0.59 

SBS 1116+583A 10.52 3.15 7.44 0.71 

Mrkl42 47.89 35.64 12.73 0.27 

NGC5548 20.86 12.42 8.37 0.40 



Note. — Col. (1): object name. Col. (2): the total continuum luminosity at 
5100 A. Col. (3): the AGN luminosity estimated from the power-law continuum 
fit. Col. (4): the host-galaxy luminosity estimated from the host-galaxy template 
fit. Col. (5): the host-galaxy fraction. 



To obtain the AGN continuum luminosity (nuclear luminos- 
ity, Lsioo.n), the host-galaxy contribution to the total luminos- 
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Fig. 2. — The host-galaxy luminosity at 5100 A, measured from each SE 
spectrum based on the spectral decomposition analysis, as a function of S/N. 
A few low-S/N spectra (open circles) are removed from further analysis to 
avoid biases due to low-quality data. Measurement uncertainties estimated 
from the method given in §3.2.3 are expressed by vertical error bars. 



been subtracted. 

Figure|2]shows the starlight luminosity measured from each 
SE spectrum as a function of S/N. As expected, the starlight 
is not constant due to the effects of seeing and miscentering. 
The variability ranges from 10% to 20% with an average of 
0.06 ±0.01 dex. These results underscore the importance of 
subtracting the host-galaxy starlight in making the rms spec- 
tra. Otherwise, the rms spectra may contain a contribution 
from the variable amount of starlight observed through the 
slit (see § 3.3 and Figure|5]l. 

As a consistency check, we directly compare the host- 
galaxy flux of NGC 5548 measured from our spectral decom- 
position with that from the HST imaging analysis as similarly 
done by Bentz et al. (2009b). In order to calculate the amount 
of light observed through the spectroscopic aperture, we used 
an aperture size of 4" x 10"1 as used in the LAMP spec- 
troscopy analysis, after smearing the point-spread-function 
(PSF) subtracted HST image with a 2" Gaussian seeing disk. 
The host-galaxy flux of NGC 5548 based on the spectral de- 
composition is 2.47 X 10"'^ erg s"' cm"^ A"', while the HST 
imaging-based galaxy flux is 2.73 x 10"'^ erg s"' cm"^ A~'. 
Thus, the small difference (~ 10%) between the two analy- 
ses shows the consistency in host-galaxy flux measurements. 
When we use a smaller seeing disk (e.g., a 1"5 Gaussian 
disk), the host-galaxy flux measured from the HST imaging 
analysis increases by ~ 13%, indicating that the actual seeing 
size will slightly change the host-galaxy flux measurements. 



\ \ \ ^ \ \ \ ^ \ \ \ 

NGC 4748 unweighted 
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4700 4800 4900 5000 

Rest Wavelength (A) 

Fig . 3 . — Comparison of the rms spectra of NGC 4748 generated with three 
different methods: unweighted rms (black), S/N weighted (red), and maxi- 
mum likelihood (blue). For this object, two bad epochs with low-S/N data 
were removed as described in § 3.2. 

ity should be subtracted from the measured total luminosity. 
In principle, the host-galaxy luminosity can be determined by 
separating a stellar component from a point source using sur- 
face brightness fitting analysis based on a high-resolution im- 
age. Such an analysis is in progress based on the HST WFC3 
images of the LAMP sample (GO-11662, PI. Bentz). For 
this paper, however, we used the information obtained from 
the spectral decomposition. We note that, although the host- 
galaxy flux should be constant, the amount of host-galaxy 
contribution to the total flux can vary in each epoch's spec- 
trum because of seeing variations and miscentering in the slit. 
Thus, the nuclear luminosity, Lsioo.n, needs to be estimated for 
each individual spectrum from which Fe II and starlight have 



3.2.3. Error Estimation 

To estimate the uncertainties of the line-width and lumi- 
nosity measurements from SE spectra, we adopted the Monte 
Carlo flux randomization method (e.g., Bentz et al. 2009(5 
IShen et al.H20Tli) . First, we generated 50 mock spectra for 
each observed spectrum by adding Gaussian random noise 
based on the flux errors at each spectral pixel. Then we mea- 
sured the line widths and AGN luminosities from the simu- 
lated spectra using the method described in §3.2. 1 and §3.2.2. 
We adopted the standard deviation of the distribution of mea- 
surements from 50 mock spectra as the measurement uncer- 
tainty. For a consistency check, we increased the number of 
mock spectra up to 100 and found that the results remain the 
same. In the case of the total luminosity (Lsioo.t), we mea- 
sured the uncertainty as the square root of the quadratic sum 
of the standard deviation of fluxes and average flux errors in 
the continuum-flux window. 

3.3. Mean and RMS Spectra 

In this section, we describe the process of generating mean 
and rms spectra, and present the method for measuring the 
line width and continuum luminosity. The mean spectra are 
representative of all single spectra, thus they are useful to con- 
strain the random errors of measurements from single-epoch 
spectra. In contrast, reverberation mapping studies generally 
use rms spectra to map the geometry and kinematics of the 
same gas that responds to the continuum variation. By com- 
paring the line profiles between rms and single spectra, one 
can investigate any systematic differences of the correspond- 
ing line widths, and therefore improve the calibration of BH 
mass estimators. 

3.3.1. Method 
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We generated mean and rms spectra for each object using 
the following equations: 



1 ^ 



(1) 



rms (A) = 



1 

i=l 



where //(A) is the flux of i-th SE spectrum (out of spectra). 

The unweighted rms spectra can be biased by low-S/N 
spectra, often showing peaky residual features in the contin- 
uum. These spurious features in continuum can affect the 
wings of the emission lines and therefore the measurement 
of line dispersion. To mitigate this effect it is best to con- 
sider more robust procedures. We considered the following 
two schemes. First, we used the S/N as a weight, with the 
following equations: 



(r(A)) = ^w,/;(A), 



(3) 



rms'^CA) = 



\ 



1 



N 

1=1 



Y.''i[fii\)-{rw)]\ (4) 



2 ;=i 



(5) 



where w,- is the normalized S/N weight defined by 
^'■ = 1^ ■ 

!=1 

Alternatively, we considered the maximum likelihood 
method. Assuming Gaussian errors, the logarithm of the like- 
lihood function (up to a normalization constant) is given by 

-2^1netot,/W-2^-^ ^^T^^, (6) 

!=1 (=1 



2lnC-- 



t,,(A) 



where 



(7) 



and e,(A) is the error in the flux fi{X). Here (/(A)) is the mean 
flux while rms(A) is the intrinsic scatter - that is, the rms flux 
after removing measurement errors. By maximizing the log- 
likelihood, we obtain the mean and rms spectra. The max- 
imum likelihood method also provides proper errors in the 
rms spectra. We calculated errors in the inferred mean and 
rms spectra in a standard way, by computing their posterior 
probability distribution after marginalizing over the other pa- 
rameters. We adopt I- a errors as symmetric intervals around 
the posterior peak containing 68.3% of the posterior probabil- 
ity. 

Figure[3]compares rms spectra of NGC 4748 generated with 
the unweighted rms method, the S/N weighted method, and 
the maximum likelihood method, after removing two bad 
epochs as described in § 3.2. As expected, the S/N weighted 
rms spectrum is less noisy than the unweighted rms spectrum. 
The rms spectrum based on the maximum likelihood method 
is similar to but slightly noisier than the S/N weighted spec- 
trum. In particular, the maximum likelihood method gener- 
ates noisy patterns around the [O III] line region, presumably 



due to the fact that the error statistics have changed owing to 
the subtraction of the strong [O III] line signals. In the case 
of the mean spectrum, all three methods produce almost iden- 
tical results. Thus, we choose the S/N weighting scheme to 
generate the mean and the rms spectra, and adopt the errors 
of the rms spectra from the maximum likelihood method. We 
note that using the rms spectra based on the maximum likeli- 
hood method does not significantly change the results in the 
following analysis. If more bad epochs with low S/N are re- 
moved in generating rms spectra (as practiced in the reverber- 
ation studies; e.g., Bentz et al. 2009c), the difference among 
the three methods tends to be smaller 

We note that there may be a potential bias in the S/N 
weighted method owing to the fact that in the high contin- 
uum state the S/N is higher while emission lines are narrower 
Thus, the S/N weighted rms spectra can be slightly biased to- 
ward having narrower lines. On the other hand, the time lag 
between the luminosity change and the corresponding veloc- 
ity change will reduce the bias since the high luminosity and 
the corresponding narrow line width are not observed at the 
same epoch. 

To test this potential bias, we compared the line-width mea- 
surements based on S/N weighted and unweighted rms spec- 
tra, respectively. We find that the line width decreases by 
2.6 ± 2. 1 % for anp and 2.7 ± 1 .3% for FWHMh/3 when the 
S/N weighted rms spectra are used, indicating that the bias is 
not significant for our sample AGNs. However, this offset is 
not due to the luminosity bias since the S/N ratio does not cor- 
relate with AGN luminosity. Instead, the change of the S/N 
ratio is mostly due to the effects of seeing and miscentering 
within the slit since different amounts of stellar light were ob- 
served within the slit on different nights. Considering the low 
level of luminosity variability and the time lag, the night-to- 
night seeing and weather variations would be the predominant 
factors affecting the S/N ratio. The average offset of ^^3% is 
dominated by two objects, NGC 6814 (0.07 dex for anp, 0.02 
dex for FWHMh/j) and SBS 1 1 16 (0.04 dex for aufs, 0.05 dex 
for FWHMh;9), which showed the largest stellar fraction in 
Fig. 1, thus supporting our conclusion. By excluding these 
two objects, the average offset decreases to ^1%. Thus, we 
conclude that the potential AGN luminosity bias in the S/N 
weighted method is not significant, at least for our sample. 

3.3.2. The Effect of Host Galaxy, Fe II, and He II 

Although the rms spectra are supposed to contain only vary- 
ing components of AGN spectra, residuals of narrow lines 
(e.g., [O III]) are often present due to residual systematic er- 
rors (due to calibration issues; Bentz et al. 2009c). Addition- 
ally, the variation of the host-galaxy starlight contribution to 
the total flux can be significant (10-20%) in the extracted SE 
spectra as discussed in §3.2.2. This variable starlight is re- 
sponsible for the stellar absorption features often visible in 
the rms spectra. To demonstrate the presence of stellar ab- 
sorption lines in the rms spectra, we fit the continuum with a 
stellar-population model. As shown in Figure |5] it is clear 
that the rms spectra show stellar absorption lines [such as 
the Mg b triplet (-5175 A), Fe (5270 A), and possibly 
(4861 A)] for Seyfert 1 galaxies having strong starlight con- 
tribution. Thus, for AGNs with high starlight fraction, like 
the ones considered here, it is important to remove the vari- 
able starlight in order to generate pure AGN rms spectra and 
correctly measure the widths of the broad lines. 

To minimize these residual features in the rms spectra, we 
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Fig. 4. — Left: S/N weighted mean spectra of 9 Seyfert galaxies. Right: S/N weiglited rms spectra. In each panel, red lines represent spectra obtained after 
removing narrow lines, Fe II emission. He II lines, and host-galaxy starlight from each individual SE spectrum. Black lines represent rms spectra obtained without 
removing the same components from each individual spectrum. Shaded regions show the eiTors from the maximum likelihood method described in §3.3. 
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Fig. 5. — The rms spectra of 4 Seyfert 1 galaxies with strong stellar fea- 
tures. Black (red) solid lines represent rms spectra (stellar model fit). Strong 
absorption lines are labelled with arrows. 

subtracted the narrow lines in all SE spectra before making 
mean and rms spectra. We also subtracted the Fe II emission 
blends, host-galaxy starlight, and the He II emission line from 
each SE spectrum. In Figure [H we show the S/N weighted 
mean and rms spectra with and without prior removal of 
narrow-line components, the Fe II blend. He II lines, and host- 
galaxy starlight. Clearly, the rms spectra are significantly af- 
fected by this procedure. In particular, removing the Fe II and 
He II emission changes the continuum shape around Hf3. For 
the objects with higher starlight fraction, stellar Hf3 absorp- 
tion is present in the rms spectra, if starlight is not removed 
from each SE spectrum. 

To quantify the change of the line widths due to prior re- 
moval of the starlight. He II, and Fe II components, we com- 
pared the line-width measurements from rms spectra gener- 
ated with/without prior removal. Each panel in Fig. |6] shows 
the effects of individual components by comparing the line- 
width measurements from the rms spectra with prior removal 
of all three components (i.e., host-galaxy stellar features. 
He II, and the Fe II blend) with those from the rms spec- 
tra without subtracting one of the three components, respec- 
tively. We found that the effect of host-galaxy stellar features 
is stronger than those of He II and the Fe II blend. Without 
subtracting host-galaxy stellar features, the rms line widths 
decrease by 18 ± 5% for cth/s and 4 ± 3% for FWHMh^, indi- 
cating that the line wings are more affected than the line core. 
The large increase of aup can be understood as the buried 



H/3 line wings within the residual of stellar features are re- 
stored by subtracting the host stellar lines, leading to a lower 
continuum level and larger line width. 

The subtraction of He II changes the line width of objects 
that show strong blending with the H/3 line (e.g., Mrk 1310 
and NGC 6814). On average, the effect of He II on the H/3 
Hne width is at the 4.7 ±2.2% level for anp and the 1 .6 ±0.6% 
level for FWHMh^. In the case of the Fe II subtraction, the 
effect on the rms line widths is more complex. Line widths 
increase for some objects and decrease for other objects, de- 
pending on whether the Fe II emission residual is strong. For 
example, if the Fe II residual is prominent in the continuum re- 
gion (i.e., 5080-5550A), then the removal of Fe II will lower 
the continuum level, increasing the H/3 line width. In con- 
trast, if the Fe II residual is strong under H/3, then the H/3 line 
width will decrease by subtracting Fe II. On average, the ef- 
fect of Fe II on the HjS Hne width is at the 1 .2 ± 2.9% level 
for (Tup and the 1 .6 ± 2.2% level for FWHMh^. 

Without prior removal of all three components (i.e., host- 
galaxy stellar features. He II, and the Fe II blend), the line 
widths are underestimated by 1 8 ± 6% for crup and 5 ± 4% for 
FWHMh^, due to the combined effects as described above. 
Subtracting stellar features has the most significant impact on 
the measurements of rms line dispersion, demonstrating the 
importance of prior removal of starlight when stellar contribu- 
tion is significant. Moreover, in order to successfully remove 
the He II blending in the rms spectra, the host-galaxy compo- 
nent as well as Fe II emission blends should be simultaneously 
fitted in the modeling of the continuum. Thus, we conclude 
that for AGNs with strong host galaxy starlight, strong Fe II, 
or blended He II, it is necessary to remove all non-broad-line 
components from SE spectra in order to generate the cleanest 
rms spectra and reduce errors in measuring the H/3 line width. 

3.3.3. Mean Spectra 

We generated the S/N weighted mean spectra without prior 
removal of narrow lines, iron emission, and host-galaxy 
starlight. Then, we used the same multi-component spectral 
fitting procedure as used for the SE spectra (see Fig.[T]). Note 
that in the case of mean spectra, removing the narrow lines, 
Fe II blends, and host-galaxy absorption features before or af- 
ter generating the mean spectra results in almost identical H/? 
broad-line profiles. 

3.3.4. Error Estimation 

Using the S/N weighted rms and mean spectra, we mea- 
sured the widths of the H/3 line from the continuum- 
subtracted spectra and determined the continuum luminosity 
at 5100 A, as described in §3.2.1 and §3.2.2. We estimated 
the uncertainty of the line-width measurements in the S/N 
weighted mean and rms s pectra using the bootstrap method 
(e.g., P eterson et al.ll2004l) . One thousand samples per object 
were generated. The median and standard deviation of the 
distribution of measurements were adopted as our line-width 
estimate and uncertainty; these are listed in Table [3] We also 
estimated the line-width uncertainties for the rms spectra us- 
ing the method given in §3.2.3. We found that the errors esti- 
mated from both the Monte Carlo flux randomization and the 
bootstrapping were consistent within a few percent on aver- 
age, which yielded almost identical fitting results. 

4. ANALYSIS AND RESULTS 

4. 1 . Testing the Assumptions ofSE BH Mass Estimators 
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Fig. 6. — The effect of each blended component on the rms line widths us- 
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without prior removal of each component (given in the upper-left corner) are 
plotted as a function of final rms line widths (all non-broad-line components 
removed). The dashed line indicates an identity relationship. The average 
offsets with their standard eiTors are given in the lower-right comer of each 
panel. 



Single-epoch Mbh estimates are based on the "virial" as- 
sumption and on the empirical relation between BLR size and 
AGN luminosity. Since Mbh does not vary over the time scale 
of our campaign, AGN luminosity and line velocity should 
obey the relation oc L"° In this section, we test this as- 
sumption by studying the relation between the line width and 
continuum luminosity from individual SE spectra of Arp 151, 
the object with the highest variability during the LAMP cam- 
paign. 

In Figure |2] we present the time variation of the line width 
and luminosity of Arp 151. Line width and luminosity are 
inversely correlated, although the variability amplitude is 
smaller in luminosity than in line width. Ideally, the luminos- 



TABLE 3 

Rest-Frame Broad H/3 Line-Width Measurements 



Object 


Mean Spectrum 


RMS Spectrum 




0"H8 


FWHMh^ 


0"H/3 


FWHMhs 




(kms-1) 


(km s-l) 


(km S-l) 


(km S-' j 


Arp 151 


1726 ±17 


3076 ±39 


1295 ±37 


2458 ± 82 


NGC4748 


952 ±6 


1796 ±8 


791 ±80 


1373 ±86 


Mrkl310 


1229 ±12 


2425 ± 19 


921 ±135 


1823 ±157 


Mrk202 


1047 ±8 


1787 ± 15 


962 ±67 


1794±181 


NGC4253 


1232 ±9 


1946 ±10 


538 ±92 


986 ±251 


NGC6814 


1744 ±12 


3129 ±14 


1697 ±224 


2945 ±283 


SBS 1116+583A 


1460 ±23 


3135 ±36 


1550±310 


3202 ±1127 


Mrk 142 


970 ±5 


1671 ±6 


700 ±54 


1601 ±224 


NGC5548 


4354 ±25 


12402±111 


3900 ± 266 


12539 ±1927 



ity variability should be four times as large as the line-width 
variability (0.042 ±0.001 dex and 0.027 ±0.003 dex, respec- 
tively, for FWHMh^ and cth/s). However, one must take into 
account the residual contamination from nonvariable sources 
to the observed continuum. In fact, the amplitude of the lumi- 
nosity variability is significantly smaller than expected based 
on the line-width variability if the total luminosity is used 
(bottom panel). In contrast, providing a validation of our con- 
stant continuum subtraction procedure, the variability ampli- 
tude of the nuclear continuum is consistent with that expected 
from the line width, as we will further quantify below. 

In Figure |8] we compare measured luminosities and line 
widths in order to test whether they obey the expected relation 
oc Lr^-^ . The continuum variation is shifted by the mea- 
sured time lag, 4 days, to account for the time delay between 
the central engine and BLR variations and then matched with 
the corresponding epochs of line-width variations. Note that 
densely sampled light curves are required for this correc- 
tion. As expected, the observed correlation between total flux 
^5ioo,t and line width is steeper than the theoretical correla- 
tion. In contrast, the correlation between nuclear flux and line 
width is consistent with the theoretical expectation. The best- 
fit slopeflare -1.46 ±0.31 (with intrinsic scatter 0.05 ±0.01 
dex) for FWHMh^, and -1.09 ±0.15 (with intrinsic scatter 
0.02 ±0.01 dex) for cth/?, which is consistent with the ex- 
pected value of -1. The linear correlation coefficients be- 
tween the nuclear luminosity and the line widths are -0.86 
for the line dispersion and -0.77 for the FWHM, indicating 
the tighter inverse correlation of continuum luminosity with 
the line dispersion than with the FWHM. 

The agreement of the observed correlations with those ex- 
pected for an ideal system is remarkable, considering the 
many sources of noise in the observed velocity-luminosity re- 
lation. They include residual errors in the subtraction of the 
host-galaxy starlight contribution and the measurement un- 
certainties of line widths and luminosities. The inverse corre- 
lation between line width and luminosity further corroborates 
the use of SE mass estimates ([Pe terson & Wandell 19991 120001 : 
lKollatschnvll2003l: iPeterson eta l. 2004). 

* We used the Bayesian linear regression routine linmix_err devel- 
oped by Kelly (2007) in the NASA IDL Astronomy User's Library. This 
method is currently the most sophisticated regression technique, which takes 
into account intrinsic scatter and nondetections as well as the measurement 
errors in both axes, generating the random draws from posterior probabil- 
ity distribution of each parameter for the given data using MCMC sampling. 
In this study, we take best-fit values and uncertainties of parameters as the 
median values and ±1ct standard deviation of 10,000 random draws from 
corresponding posterior distributions. 
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Fig. 7. — Time variations of the H/3 line width (top; FWHMnij; uppermid- 

dle: (Jnfj) and the continuum luminosity at 5100 A (lower middle: nuclear; 
bottom: total) of Arp 151. The dashed lines represent the average values over 
the monitoring period. The rms dispersion values are given in each panel. 

4.2. Uncertainties Due to Variability 

Since the line width and continuum luminosity of an AGN 
vary as a function of time, mass estimates from SE spectra 
may also vary. Owing to its stochastic nature, this variability 
can be considered a source of random error in SE mass esti- 
mates. In this section, we quantify this effect by comparing 
SE measurements with measurements from the mean spectra. 

4.2.1. The Effect of Line-Width Variability 

We quantify the dispersion of the distribution of line-width 
measurements using all SE spectra. This dispersion can be in- 
terpreted as a random error due to the combined effect of vari- 
ability and measurement errors. In Figure |9] we present the 
distributions of FWHMh^ measurements from all SE spectra, 
after normalizing them by the measurement from the mean 
spectra. All SE values are normalized to the FWHM mea- 
sured from the mean spectra. The standard deviation of the 
FWHM distributions ranges from 0.009 dex to 0.042 dex, 
with an average of 0.021 ±0.004 dex (^5%) across all ob- 
jects. Note that the standard deviation includes the variability 
and the measurement error 

In Figure[TO] we plot the distributions of line dispersion for 
all objects. The dispersion of distributions ranges from 0.013 
dex to 0.040 dex, with an average and rms of 0.023 ±0.003 
dex (^^5%) for the entire sample. SBS 1116 shows the broad- 
est distribution; however, part of this scatter can be attributed 
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Fig. 8. — Test of SE mass estimates for Aip 151. Each filled circle repre- 
sents SE measurements after shifting luminosity measurements by the aver- 
age lag of 4 days (left: total luminosity; right: nuclear luminosity). Dashed 
fines represent the correlation L"^ oc V^^ expected from the virial theorem 
and the size-luminosity relation, while solid lines are the best-fit slopes. The 
values of best-fit slopes and its uncertainties are given in each panel. 
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to the residual systematic in the left wing of H/? due to the bad 
pixels in the original spectra, as discussed previously. 

By averaging the standard deviation of the distribution of 
the line-width measurements for all 9 objects in the sample, 
we find that the uncertainty of SE BH mass estimates due to 
the line-width variation and measurement errors is on aver- 
age 0.044 dex. Note that the dispersion of the line-width 
distribution strongly depends on the variability. For exam- 
ple, Arp 151 has the largest variability amplitude and also the 
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Fig . 11 . — Distrib ution of nuclear luminosities measured from single-epoch 
spectra (see j|3.2.2) . Each luminosity value is normalized to the nuclear lu- 
minosity measured from the mean spectra. The standard deviation of the 
distribution is given in each panel. The average standard deviation of all 9 
objects is 0.048 ± 0.008 dex. 



largest variability in the line width. This is expected if line 
flux correlates with BLR size and both are connected to the 
BH mass. Based on these results, we conclude that the typical 
uncertainty of SE mass estimates due to line-width variabil- 
ity is ^10%. However, as discussed below, this uncertainty is 
partly cancelled out in the virial product by the inverse corre- 
lation with the variability of the continuum. 



X (7^^ as measured from 



the nuclear luminosities at 5100 A, after normalizing them 
by the nuclear luminosity measured from the mean spectra. 
The standard deviation of the luminosity distributions ranges 
from 0.019 to 0.097 dex, with an average of 0.048 ± 0.008 
dex (^12%), which can be treated as a random error of the 
continuum luminosity measured from a SE spectrum due to 
the luminosity variability and measurement error. 

Based on the empirical size-luminosity relation, the random 
errors of the luminosity enter the uncertainty of the SE mass 
estimates as the square root (i.e., 0.024 dex). This is some- 
what smaller than the uncertainty of SE mass estimates due 
to the line-width variability, 0.044 dex, as determined in the 
previous section, indicating that the two do not cancel each 
other exactly. 

4.2.3. Combined Effect 

Since the luminosity and the line width are inversely corre- 
lated as cx L~"^, one may expect that the variability of lu- 
minosity and line width cancel out in the SE mass estimates. 
However, the two effects may not compensate each other ex- 
actly, for a variety of reasons. First, there is a time lag between 
continuum and emission-line variability. Second, variations 
such as in the ionizing flux may indicate that the luminosity 
at 5100 A traces the broad-line size only approximately. In 
order to quantify the combined effect of the continuum lumi- 
nosity and line-width variability, we thus investigate the dis- 
tribution of the virial product ^5100 « ^ 
SE spectra. 

In Figure [121 we present the distribution of the SE virial 
products, normalized by the virial product measured from the 
mean spectra. The standard deviation of the distributions can 
be treated as a random error due to the combined variabil- 
ity and measurement errors. The average rms scatter (cor- 
responding to a source of random measurement errors when 
using the SE estimator) of the virial products is 0.052 ±0.006 
dex when the line dispersion (cth/s) is used, and 0.049 ± 0.006 
dex when FWHMh^ is used. 

In agreement with previous studies jWiUiite et al.l 120071 : 
IWoo et all 120071: iDeniiev et all l2009h . fliese results suggest 
that BH masses based on SE spectra taken at different epochs 
are consistent within ^ 0.05 dex (^ 12%) uncertainty, negli- 
gible with respect to other sources of uncertainty which are 
believed to add up to ^0.4-0.5 dex (see §5.1). 

4.3. Systematic Difference between SE and Reverberation 

Masses 

In order to assess the accuracy of the SE mass estimates, we 
need to compare the SE masses with the masses determined 
from reverberation mapping. Setting aside potential differ- 
ences in the virial coefficient, there are two main sources of 
systematic uncertainties in SE mass estimates. One is the po- 
tential difference of the line profile between SE spectra and 
the rms spectra. The other is the systematic uncertainty of the 
size-luminosity relation. We postpone discussion of the lat- 
ter to a future paper when more accurate HST-hased nuclear 
luminosities will be available. Therefore, in this section we 
focus on the systematic difference of the H/3 line profile and 
derive new SE mass estimators recalibrated to account for the 
difference found. 



4.2.2. The Effect of Luminosity Variability 

We now consider the effect of luminosity variability on SE 
mass estimates. In Figure (TT] we present the distributions of 



4.3.1. Comparing Line Profiles 

In Figure [13] we compare the broad H/3 line profiles mea- 
sured from the mean and rms spectra after normalizing by 
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the peak flux. Generally the H/3 line is broader in the mean 
spectra than in the rms spectra, indicating that the varia- 
tion is weaker in the line wings than in the line core. It is 
worth noting that the observed offset cannot be explained by 
the contamination of the narrow H/3 component or the Fe II 
blends since we consistently subtracted them in both the mean 
and rms spectra. To verify this we arbitrarily decreased the 
amount of narrow component subtracted from the observed 
H/3 profile, and found that the large offsets between rms and 
mean spectra are virtually unchanged. 

The broader line width in the mean spectra has been noted 
in previous reverberation studies (e.g., Serge ev et al. 1999; 
Shapovalova et al. 2004). ICoUin et al.l (1200 6) reported that 
the line widths in the mean spectra wer e typically bro a der by 
~ 20% than those in the rms spectra. iDennev et al.l (120101) 
also found that some objects in their reverberation sample 
clearly showed narrower line widths in the rms spectra than 
in the mean spectra. Several different and somewhat mutu- 
ally exclusive explan ations have been sug gested for this dif- 
ference. For example. [Shields et al.l(ll995h explained the sys- 
tematic difference of the line width as being due to the high- 
velocity gas in the inner BLR being optically thin to the ion- 
izing continuum and hence fully ionized. In this way, the 
line wings have weak variability and are sup pressed in the 
rms spectra. In contrast, Korista & GoadI (120041) suggested a 
distance-dependent responsivity of optically thick clouds to 
explain the weak variability of Balmer line wings. 

We quantify the systematic offset in line width in Figure [14] 
by showing the ratios of the line width measured from the 
mean (and SE) spectra to those measured from the rms spectra 
as a function of line width. The average offset in FWHMh/j 
is 0.07 ±0.03 dex (0.05 ±0.02 dex, if NGC 4253, the object 
with the narrowest line, is excluded). In the case of line dis- 
persion ((Jhis), the offset is slightly larger, 0.10 ±0.04 dex 
(0.07 ±0.02 dex if NGC 4253 is excluded). The larger offset 
of the fine dispersion in comparison with FWHM is consis- 
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Fig. 13. — Comparison of the H/3 broad-line profiles in the mean (black) 
and rms (red) spectra. In each panel, solid lines represent the data, while dot- 
ted lines represent the Gauss-Hermite series fitting results. Each line profile 
is normalized by the maximum value of the fit. Dashed vertical lines indicate 
the center of the H/3 line. 



tent with there being mainly a difference between variability 
in the wings and in the core. 

There seems to be a systematic trend, in the sense that the 
offset becomes relatively larger for the narrower line objects, 
but its origin is not clear In particular, the narrow-line Seyfert 
1 galaxy NGC 4253 (Mrk 766) has the narrowest H/3 line in 
the sample and shows the largest systematic difference. It is 
possible that the systematic difference for this particular ob- 
ject with very narrow H/3 (FWHMh/3 (rms) < 1000 km s"') 
may be amplified due to imperfect subtraction of the narrow 
component, the Fe II blends, or starlight. However, the trend 
is present even if we remove this object from the sample. 

In order to correct for this potential bias, we derive a re- 
lation between H/3 line width as measured from rms and SE 
spectra by fitting the trend as shown in Figure [141 Using the 
linear regression routine linmix_err (Kelly 2007), we fit 
the linear relationship in log-scale using bootstrap errors de- 
termined in §3.3. We also determined the slope excluding the 
narrowest-line object (NGC 4253) or the broadest-line object 
(NGC 5548) from the sample. As shown in Figure 13, remov- 
ing either NGC 4253 or NGC 5548 from the sample does not 
significantly change the slope. 

In addition, we fit the slope using a fixed error for all ob- 
jects. Since the bootstrap errors on the rms line widths are 
significantly different for each object owing to the different 



Recalibrating Single-Epoch Black Hole Mass Estimates 



13 





b 



0.3 



c 

CO 



O 



0.0 



H 1 1 1 h 



I I I I 



_ i _i ^ 



3.0 



3.5 



4.0 



log FWHMh7 (km/s) 



3.0 3.5 

log cj-^ (km/s) 



Fig. 14. — Direct comparisons of line widths and their ratios, i.e., FWHMh/3 (left) and (Th^ (right), measured from SE (or mean) and rms spectra, as a function 
of line width. Average offset is 0.07 ± 0.03 dex for FWHMh/3 and 0. 10 ± 0.04 dex for cth/j . Dashed Hnes indicate an identity relation while solid lines are the 
best-fit results using all objects (red), excluding NGC 4253 (blue), or excluding NGC 5548 (brown) using bootstrap errors. 



quality and S/N ratios of individual single-epoch spectra, we 
assigned a fixed error, such as 20% on both axes, to test the ef- 
fect of errors. The best-fit slope using a fixed error is slightly 
shallower than that with bootstrap errors since the most off- 
set object, NGC 4253 has a large bootstrap error and conse- 
quently has smaller weight in the fitting process. 

To secure a large dynamic range, we decided to use all 9 
objects for the fit and adopt the best-fit result using bootstrap 
errors. The adopted best fits are expressed as 



logFWHMH/3(rms) = 

-0.405(±0.051)+ 1.095(±0.015) logFWHMH/3(SE), 



(8) 



logcrH^)('-«M) = -0.434(±0.060)+1.106(±0.019)logcrH^)(5£). 

(9) 

However, these fits should not be extrapolated to high- 
velocity objects; otherwise, negative bias will be introduced. 
Since our sample consists of relatively narrow-line Seyfert 1 
galaxies, we recommend that readers use Eqs. [8]and|9]for ob- 
jects with FWHMh/3,se < 3,000 km s"' and cth^^se < 2,000 
km s"', respectively. 

4.3.2. Systematic Ojfset of Mass Estimates 

The systematically broader line width in SE spectra would 
result in overestimates of SE masses if unaccounted for. In 



Figure [15] we plot the ratio of the SE virial product (VP) with 
respect to the virial product based on the reverberation stud- 
ies as a function of the virial product. Note that to demon- 
strate the effect of the systematic difference between SE and 
rms spectra, we simply used the measured /?blr for all SE 
virial products, instead of using L5100 and the size-luminosity 
relation. As expected, virial products exhibit a biased sys- 
tematic trend. The average offset of all SE measurements is 
0. 152 ± 0.009 dex when FWHM is used in the virial product 
(top), and 0.204 ±0.011 dex for the dH^-based virial prod- 
uct (bottom). The average offset of all measurements based 
on mean spectra is 0.147 ± 0.066 (0.201 ± 0.075) dex when 
FWHMh^ (ch^) is used in the virial product. To avoid po- 
tential biases from the narrowest-line object, we recalculate 
the average offset after removing NGC 4253. The average 
offset of all SE measurements is now 0.093 ± 0.006 dex and 
0. 135 ±0.007 dex for FWHMnfi and cth/s, respectively. When 
the measurements from mean spectra are used in comparison, 
the average offsets are 0.092 ± 0.040 dex and 0. 136 ± 0.043 
dex for FWHMh/S and aiip, respectively. Thus, the SE BH 
masses of broad-line AGNs with virial products in the range 
^10^"^ Mq can be overestimated by ^25-35% if the same 
recipe used for rms spectra is adopted. 

These results are similar to the findings by Collin et al. 
(2006), who investigated the systematic offset of vmal prod- 
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(FWHMna) is used for the virial products. 

uct estimates between rms and mean spectra using a different 
sample of reverberation mapping measurements. Although it 
is not straightforward to directly compare their results with 
ours since the methods of generating rms spectra and mea- 
suring line widths are substantially different, the similar sys- 
tematic offset between SE and reverberation masses clearly 
demonstrates the importance of calibrating SE masses. 

4.3.3. Comparing FWHMh/i and uwp 

We compare FWHMh^ and a^p in Figure 15. As previ- 
ously noticed in other studies (e.g., Collin et al. 2006; McGill 
et al. 2008), the shape of the H/3 line is different from a Gaus- 
sian profile, for which FWHMh^/cth,? is expected to be 2.35. 
As shown in Figure 16, narrower lines tend to have stronger 
wings leading to a lower FWHMh^/cth/? ratio, while broader 
lines are more core dominated with a higher FWHMh^/cth^ 
ratio. These results are consistent with those of previous stud- 
ies, although our sample is composed of objects with narrower 
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Fig. 16. — Comparisons of FWHMh/3 and cth^, using measurements from 
rms spectra (top) and SE spectra (bottom). The dashed line indicates an iden- 
tity relationship while the dotted line represents the Gaussian profile, i.e., 
FWHMna/(Tn^ = 2.35. The red solid Une shows the best-fit relation. 

lines than previously studied. 
The best-fit correlation based on the rms spectra is 

logcrH^(rms) = 0.401(±0.395)-l-0.792(±0.117) logFWHMH^j(rms). 

(10) 

In the case of measurements from the SE spectra, the best fit 
is expressed as 

logo-H/^(SE) = 0.567(±0.027)-l-0.753(±0.008) logFWHMH/^(SE). 

(11) 

We note that these results are somewhat limited by the small 
dynamic range of our sample and the lack of objects with 
FWHM > 3000 km s"'. Further analysis with broader line 
objects is required. However, in the case of broader line ob- 
jects in the literature, we do not have consistently measured 
line widths from rms spectra as described in §3.3. Neverthe- 
less, we will use this fit to convert FWHMh;3 to cth^ in §4.3.4. 

4.3.4. Line-Width Dependent Mass Estimators 

In order to avoid potential systematic biases in SE spec- 
tra, we derive line-width dependent mass estimators, using the 
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best-fit relations derived above (see Fig. [I4] |. As a reference, 
we use the mass estimator normalized for the virial product 
from rms spectra (reverberation results) using the virial factor 
log / = 0.72 determined from the Mbh - c* relation of rever- 
beration mapped AGNs (Woo et al. 2010): 



Mbh = 10^-^"^Me 



fTH,3(rms) \ / AL: 



lOOOkms"'/ V 10^4 erg s- 



0.518 



(12) 



If we replace (Th/s from rms spectra with aap measured from 
SE spectra using Eq. |9] the mass estimator changes to 



Mbh = 10^"°M 



/ gH/3(SE) 

V 1000 km s"' 



2.212 



'5100,;i 



V 1044 erg s- 



0.518 



(13) 

As in the case of Eq. |9] we recommend readers use Eq. [T3]for 
AGNs with anp< 2,000 km s"'. 

In the case of FWHMh^, the virial factor has not been de- 
termined by Woo et al. (2010), but we can use the relations 
found above to derive a consistent expression. If we replace 
(Tup with FWHMh^ using Eq. [TOl then Eq. [TSlbecomes 

^ FWHMHMrms) y"^Y ^' 



Mbh = 10' 



V 1000 km s" 



AL- 



5100.« 



V 1044 erg s- 



(14) 

In order to use FWHMh/j measured from SE spectra, 
FWHMh;3 from rms spectra in Eq. [14] can be replaced by 
FWHMh;3 from SE spectra, using Eq. |8] Then, the mass esti- 
mator becomes 



Mbh = 10''''^''M, 



/ FWHMhMSE) \ ' "^ / AL5ioo.„ 
®V 1000 km S-' J UO^'ergs- 



0.518 



(15) 



Alternatively, we can use Eq. [13] and replace aap with 
FWHMh;3 measured from SE spectra using Eq. [TT] Then 
Eq. [T3]becomes 



^ _in6.985^, / fwhMhM^£) V'"7 ^ 

""^"-'^ 1000 km S-' ) 1 1044 ergs- J 

(16) 

which is almost identical to Eq. [15] For a consistency check, 
we compared the SE masses estimated from Eq. [T5]with those 
from Eq. [16] They are consistent within ^1%, indicating that 
Eq. 15 and 16 are essentially equivalent. As in the case of Eq. 
[8] for AGNs with FWHMh^< 3,000 km s"' we recommend 
readers use Eq. [T6]instead of Eq. [15] since the SE masses de- 
rived from Eq. 16 are slightly more consistent with the masses 
determined from Eqs. 12 and 14. 

The BH masses derived from the new mass estimators are 
consistent with each other within a ^2% offset, indicating that 
the systematic difference in the line widths between SE and 
rms spectra is well calibrated. In contrast, the ^0.2 dex scat- 
ter between various mass estimators reflects a lower limit to 
the uncertainties of our line-width dependent calibrations. In 
a sense, these new estimators can be thought of as introduc- 
ing a line-width dependent virial factor to correct for the sys- 
tematic difference of the geometry and kinematics of the gas 
contributing to the SE line profile and that contributing to the 
rms spectra. Regardless of the physical interpretation, these 
new recipes ensure that mass estimates from SE spectra and 
rms spectra can be properly compared. 

5. DISCUSSION AND CONCLUSIONS 



0.518 



5.1. Random Uncertainty 

We investigated the precision and accuracy of BH mass esti- 
mates based on SE spectra, using the homogeneous and high- 
quality spectroscopic monitoring of 9 local Seyfert 1 galaxies 
obtained as part of the LAMP project. We find that the un- 
certainty of SE mass estimates due to the AGN variability is 
^ 0.05 dex (~12%) . Our result is slightly less than that of 
iDenney et al.l ( |20 09*). who reported ^0.1 dex random error 
due to the variability based on the investigation of Seyfert 1 
galaxy NGC5548 using data covering '--^10 years. For higher 
luminosity AGNs, the uncertainty due to variability can be 
smaller since the amplitude of variability inversely correlates 
with the luminosity (e.g., Cristiani et al. 1997). For example, 
by comparing SE spectra with mean spectra averaged over 
'--^10 multi-epoch data of 8 moderate-luminosity AGN, Woo 
et al. (2007) reported that intrinsic FWHM variation of the 
H/3 line is ^1%, resulting in ~15% random error in mass es- 
timates. 

In addition to the uncertainty related to variability, the total 
random uncertainty of SE mass estimators includes the uncer- 
tainty in the virial factor, and the scatter of the size-luminosity 
relation. The scatter of the AGN Mbh - relation (Woo et al. 
2010) provides an upper limit to the random object-to-object 
scatter in the virial factor of 0.43 dex. By adding 0. 1 dex due 
to variability and 0.13 dex scatter from the size-luminosity 
relation in quadrature (assuming they are uncorrelated), the 
upper limit of the overall uncertainty of SE mass estimates is 
found to be 0.46 dex. This is consistent with the uncertainty 
of 0.4-0.5 dex estimated by Vestergaard & Peterson (2006). If 
we assume more realistically that 0.3 dex of the scatter in the 
Mbh-<^* relation measured by Woo et al. (2010) is intrinsic 
scatter (e.g., Giiltekin et al. 2010) and not due to uncertain- 
ties in the virial coefficient, then the uncertainty of the virial 
factor becomes 0.31 dex, resulting in an overall uncertainty of 
^0.35 dex in SE mass estimates. More direct measurements 
of the virial coefficient (e.g., Davies et al. 2006; Onken et al. 
2007; Hicks & Malkan 2008; Brewer et al. 201 1) are needed 
to break this degeneracy. 

Note that measurement errors in the line width and contin- 
uum luminosity are negligible in our study owing to the high 
quality of the data. However, often such high-quality data 
are not available, and measurement errors of the line width in 
particular can be a significant contribution to the total error 
budget. For example. Woo et al. (2007) estimated the propa- 
gated uncertainty in the SE mass estimates due to the FWHM 
measurement errors as 0.11 dex (30%) based on spectra with 
a S/N of 10-15. Therefore, the estimated overall uncertainty 
of ^0.35 dex should be taken as a lower limit for typical SE 
mass estimates based on optical spectra such as those from 
the Sloan Digital Sky Survey. 

5.2. Difference in Line Profile between SE and RMS spectra 

We confirmed that the H/3 line width measured from a mean 
or SE spectrum is systematically larger than that measured 
from an rms spectrum. The systematic difference corresponds 
to an average difference in virial product of ^0.15-0.20 dex. 
However, the average difference is dominated by the narrow- 
est line objects in the sample, with a decreasing trend to- 
ward broader line objects as shown in Fig. [15] (cf., Collin 
et al. 2006). These results indicate that for narrow-line AGNs 
(FWHMh/3 < 3000 km s"'), BH masses based on SE spec- 
tra can be overestimated by '^25-35% if standard recipes are 
used. In order to correct for the systematic difference of the 
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line profile, we derive new empirically calibrated line-width 
dependent SE mass estimators. 

It is important to notice that line-width measurements from 
rms spectra can also be systematically biased by residuals 
of narrow-line components, Fe II, and host-galaxy starlight. 
Fluctuations of these components can generate signatures in 
the rms spectra, resulting in a biased line profile and an im- 
proper continuum fit. We have demonstrated several new 
strategies to mitigate these effects. First, we adopt two ro- 
bust methods to derive rms spectra — using S/N weights 
and adopting a maximum likelihood approach — that mini- 
mize the contamination of low-S/N spectra when there is a 
large range in S/N. Second, we subtract Fe II and host-galaxy 
starlight based on spectral decomposition analysis of each 
individual-epoch spectrum before making the rms and mean 
spectra. These new methods substantially improve the qual- 
ity of rms spectra in measuring the width of the H/3 line for 
AGNs with strong starlight, Fe II, or blended He II. Especially 
in these cases we recommend the new methods as an useful al- 
ternative to the traditional simple methods in future reverbera- 
tion mapping studies (and possibly to revisit previous studies). 



The new methods introduced here differ from previous 
ones. The overall goal is to account for systematic uncer- 
tainties and correct for biases stemming from known effects 
such as stellar line contamination. However, it is possible 
that they might introduce biases due to unknown systemat- 
ics, especially when directly compared to other measurements 
obtained with previous techniques. An absolute comparison 
would require a third way to measure the same quantities (e.g., 
BH mass); however, we can estimate any differential bias by 
comparing our measurements to those given by Bentz et al] 
(2009c) using traditional methods. On average, (jH^i measured 
with our new methods increases by 7 ± 6% and FWHMh;? in- 
creases by 9 ±4% (see Fig. 17), suggesting that the system- 
atic uncertainties due to the new methods is lower than 10%, 
although we assume the difference is entirely caused by the 
systematic uncertainty of the new scheme. 

The average small difference between our new measure- 
ments and those of Bentz et al. (2009c) is due to the fact that 
various effects, e.g., S/N weighting, prior subtraction of host 
galaxy starlight and blended emission lines, line profile fit- 
ting, and removal of narrow H/3, are mixed together and can- 



Recalibrating Single-Epoch Black Hole Mass Estimates 



17 



celed out for individual objects. For example, if we separate 
the effect of host galaxy starlight, the difference of the line 
dispersion measurements with/without removal of host galaxy 
starlight is much la rger than the average difference between 
iBentzet al.l (I2009cl) and ours (see §3.3.2). The 10-20% rms 
scatter between previous (Bentz et al. 2009c) and new mea- 
surements indicates that the contribution of each effect varies 
on an object-by-object basis (see Fig. 17). The new meth- 
ods are useful for reducing the uncertianites of individual BH 
masses due to those various effects, and better constraining 
the intrinsic scatter of the Mbh-c* relation (e.g.. Woo et al. 
2010). When other systematic uncertainites, i.e., the virial 
coefficient, can be constrained and reduced in the future, the 
new methods will become more important for BH mass deter- 
mination. 

From the point of view of interpretation, the systematically 
narrower line width in the rms spe ctra can b e explained by 
the photoionization calculations of lKorista"& Goad (2004), 
which predict that high-velocity gas in the inner BLR has 
lower responsivity, leading to lower variability of the line 
wings and therefore a narrower profile in the rms spectrum. 
However, it is not clear why the effect should be stronger for 
the narrower line AGNs. Further investigations are required 
to reveal the nature of this systematic trend. We conclude by 
noting that our study is limited to relatively low luminosity 
and narrow line width (and hence small BH mass) AGNs. In 
future work we plan to expand our study to cover a larger dy- 
namic range in luminosity and line width. 

5.3. Implications for the Evolution of BH Host-Galaxy 
Scaling Relations 

Virtually all observational studies of the evolution of the 
BH host-galaxy scaling relations over cosmic time are based 
on SE mass estimates. Evolutionary trends are generally 
established by comparing measurement of distant samples 
(based on SE BH mass estimates) with the local scaling re- 
lations. The latter can be determined either based on SE BH 



mass estimates, or on reverberation BH mass estimates, or on 
spatially resolved kinematics. In the case of comparison be- 
tween SE mass estimates at high redshift and mass estimates 
of local objects based on different methods, it is essential that 
the mass estimates be properly calibrated. For example, the 
positive bias of H/3 line width of SE spectra compared to the 
rms spectra could lead to overestimation of BH masses for 
distant samples, if compared with local samples based on rms 
spectra. 

However, as shown in Figures 13 and 14, this bias is 
only significant for the narrower line objects with BH mass 
< 10^ Mf7^). Ty pically, high-redshift studies such as those 
by iBennert et ail dMlOj and Merloni et aL (.2010 ) focus on 
higher mass BHs, where the bias is believed to be negligi- 
ble. We note that one could completely eliminate this bias 
by comparing distant and local BH mass estimates based en- 
tirely on se lf-consistent SE BH mass estimates dWoo et alJ 
|2008HBennert et al. 201 la b). Even then, of course, one must 
keep in mind the differential nature of the measurement. For 
example, the slope infer red for the lo cal scaling relations 
based on SE spectra (e.g.. lGreene & Ho| [2006) may be biased 
with respect to the true slope if the mass estimator is biased at 
low masses, and yet one may still infer the correct evolution 
even for low masses if the bias does not change with redshift. 
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